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Abstract 

We present the exact solution for a set of nonlinear algebraic equations 
— = Trd + —JZ™^i — - — ■ These were encountered by us in a recent study 
of the low energy spectrum of the Heisenberg ferromagnetic chain . These 
equations are low d (density) "degenerations" of more complicated transcen- 
dental equation of Bethe's Ansatz for a ferromagnet, but are interesting in 
themselves. They generalize, through a single parameter, the equations of 
Stieltjes, xi = J2 m ^i V( x z ~~ x m ), familiar from Random Matrix theory. It 
is shown that the solutions of these set of equations is given by the zeros 
of generalized associated Laguerre polynomials. These zeros are interesting, 
since they provide one of the few known cases where the location is along a 
nontrivial curve in the complex plane that is determined in this work. Using a 
"Green's function" and a saddle point technique we determine the asymptotic 
distribution of zeros. 



I. INTRODUCTION 



The study of nonlinear algebraic sets of equations that arise in various physical contexts is 
a rich field of study. A famous example is the problem of Stieltjes, namely the set of equations 
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x i — J2m^i —Xm), for n variables xi, that arises when we consider the stationary points 
of the Random Matrix Gaussian Ensemble "action" 1/2J2%J — 1/2 J2i^ m ^°s(\ x i ~ x m\) [0- 
In that context, the well known result is that the x\ are all real, and further are the roots of 
Hermite's polynomials of degree n, forming a dense set along the real line with the familiar 
semicircular density of states as n — > oo. 

In a recent study of the "most elementary excitations" of the Heisenberg ferromagnet in 
1-d we came across a one parameter generalization of the Stieljes problem 

1 , 2d 1 

_ = 7T d + — }^— . (1) 

z l n m£l Z[ Zm 

where the single parameter "d" ( < d < 1) has the physical significance of the density 
of (hardcore) particles on a lattice. By a simple redefinition of variables (see below), this 
is recognizable as a generalized Stieltjes problem 1+ ^ = ~ J2 m ^i ; where, as the 
parameter A vanishes, it reduces to the standard problem. 

We found the equation quite fascinating in its own terms, requiring some new tricks to 
solve. It further appears that this problem connects with P,|8HTll that of the asymptotic 
distribution of zeros of the standard polynomials (orthogonal in most instances though not 
in the present case). It provides an explicit example of a case where the zeros live on a 
nontrivial curve in the complex plane, which is determined here. 

The paper is organized as follows. In Sec. II we obtain a perturbative solution of Eq. (|]) 
while in Sec. Ill we give the exact solution. The construction of the "Green's function" and 
the determination of distribution of zeros are described in Sec. IV. We summarize our results 
in Sec. V. Some results of Stieltjes on the properties of zeros of orthogonal polynomials are 
described in Appendix (A) . In Appendix (B) we describe a second method of obtaining the 
Green's function. 



II. THE PERTURBATIVE SOLUTION 



We first present a "perturbative solution" in a parameter displayed below, perturbing 
around the Stieltjes example. In Eq. ([[]) we make a change of variables zi = (1 + i\j3i)/(7id), 
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with A = y2d. The new variable (3i satisfies the following equation: 



1 + iXPi n^ft- f3„ 

We note that for A = 0, this equation reduces to the well known form whose solution can 
be found in terms of the roots of the Hermite polynomials, the example of Stieltjes ||. The 
general formalism for this is summarized and expanded upon in App. Thus one obtains 
Pi = xi where xi satisfy H n (y/nxi) = 0, where H n is the Hermite polynomial of degree 
n. The division by n in Eq(§) is our preference, since it enables an easier passage to the 
"thermodynamic limit"" n — > oo, and we discuss this feature later. 

One possible approach is to try for a perturbative solution around the known result of 
Stieltjes to various orders in the small parameter A (see below). The perturbation analysis 
requires a knowledge of the eigenvectors and eigenvalues of an interesting matrix studied 
by Calogero (we call this the Calogero matrix) in which the matrix elements are algebraic 
functions of xi, the roots of the Hermite polynomials. 

For A 7^ we try a perturbation series solution of the form /?/ = J2kfli (zA) fc . Putting 

this in Eq. (|1]) and matching terms order by order we get the following set of equations for 

(k) 

the unknown coefficients : 
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As noted above, the lowest order solution is given by (3^ = x\ where H n {^/nxi) = 0. It 
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turns out that it is possible to solve the equations of the perturbative series at every order. 
To see this we first define the Calogero matrix |J: 

We then note that the set of equations in Eq. (|3|), for k > 0, have the following general 
structure 

A (fc) + 1 E^4 fe) = ^ ) (5) 

n m 

(k) 

where the function g\ appearing at fcth order in perturbation is a function of lower order 
terms only and so is known. Hence at every stage we have to essentially solve a linear matrix 
equation with the same matrix appearing at all orders. This can be done by using some 
rather special properties of the Calogero matrix, T. The following result is easily obtained. 

Tmtfn = rx\ - (r ~ 1) ^~ r) x[- 2 - ^(r - 2, - 1)(± x^x^ 2 (6) 

Z s=l m=l 

where [p] denotes the largest integer < p. Note that the rhs of Eq. (^J) involves powers of xi 
alone. It is also possible to show that at every order the function gjy is a known polynomial 
of degree (k + 1) in the variable x\. Thus it follows that the fcth order solution can be 
obtained in the form 

Pl k) = J24 k H, (7) 

r=0 

where the coefficients are obtained by inserting the above form into Eq. (|5[) , using the 
property Eq. (|6|) of the Calogero matrix T and then comparing lhs and rhs of the equation. 
Till second order we get: 

Pi = xi + i\[\xf + 1(1 - I)] + (*A) 2 [^ + ^(14 - + 0(A 3 ) (8) 

It is straightforward, though tedious, to carry the perturbation to any order. We have been 
unable to find a closed form for the general term. 

As an example of this perturbation theory we note that the first order terms in A is 
already enough to compute several objects of interest. In the ferromagnet problem, one 
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needs the "energy" of the state w , defined by w = 4/ (nd) J2i l/ z i-( The present w = ne/dJ 
or Ref( []]])). To obtain this energy to order d 2 the perturbation series till first order in A 
suffices. In the limit of large n the xi form a continuum stretching from —\/2 to \/2 with 
the familiar semicircular density of states p(x) = ^\/2 — x 2 . This solution can be used to 
obtain the energy to order d 2 . The energy w for low d can be found as w = Air 2 [d — d 2 {6 < 
($ (0) ) 2 > +4 < $ 1] >} + 0(d 3 )}, where the averages are normalized sums over the indicated 
variables. Using the explicit expression Eq (H) and converting the sums to integrals over the 
semicircular density of states we get finally the low density formula 

w = A7T 2 d{l-d) +0{d 3 ). (9) 



III. EXACT SOLUTION 

The exact solution of Eq. ([[]) can be obtained using the results of Stieltjes which are 
described in App. (JAJ). We will give two different but related solutions. The first one is 
directly related to the perturbative approach. We first note that after change of variables 
Ui = y/nf3i, Eq. (H) becomes of the form Eq. ( |A8| ) with p(u) = (l + i(2d/n) 1 ^ 2 u); q{u) = —2u. 
The {ui} are then obtained as zeros of the the corresponding polynomial function g(u) which 
in this case satisfies the following differential equation: 

(1 + ipu)g"{u) - 2ug{u) + 2ng(u) = 0, (10) 



where v = ^2d/n. For v = this is just Hermite's equation and we recover our zeroth order 
perturbative result. For v ^ 0, we obtain the following series solutions (which we denote by 
Qn{u)): 

n 

Qn{u) = C kU k With 
k=0 

-(Jfe + 2)(Jfe + l) ivk(k + l) 

° k ~ 2{^k) ° k+2 ~ 2(n - k) ° k+1 

c ra - 1= -^ ( 7 1)C " ; c n = 2^. (11) 
5 



To obtain our second solution to Eq. ([!]), we make the change of variables z\ = —yi/(iin), 
so that Eq. ([[]) is transformed to: 

n/d + y l=2 ^^_ (12) 

Vi mz £i yi Um 

This is in a form where we can once again apply Stieltjes result. The corresponding dif- 
ferential equation in this case is the associated Laguerre equation Eq. ( |A11| ). Thus {z{\ 
are obtained as zeros of the associated Laguerre polynomials L(~ n l d ~ 1 \—rmz) ||. Note 
that the usual orthogonal Laguerre polynomials L^iy) have a > —1, while in our case 
a = —(n/d+1) < — (n + 1). In this case it can be proved || that for n even there are no real 
zeros while for odd n there is a single real zero. Physically this corresponds to the fact that 
the single-particle potential [see App. (^)] in this case is no longer confining and therefore 
we cannot get any position of equilibrium for the particles. 
The two solutions described are related as: 

Qn(x) = 1 ^ /d)n/2 L- n/d -\-i(2n/d) l / 2 x - n/d). (13) 

We note that since in the limit d — > we get Q n (x) — > H n (x), this leads to the following 
interesting identity relating the Hermite and Laguerre polynomials: 

lim 1 {7 ; /d)n/2 L^ d -\-i{2n/df/ 2 x - n/d) = H n (x). (14) 
Such a relation is of course well known, e.g. [0 for the case of large positive order m in 

IV. ASYMPTOTIC DISTRIBUTION OF ZEROS 

In the previous section it was shown that the solutions to Eq. (|l]) are the zeros of gener- 
alized associated Laguerre polynomials. Unlike in the case of usual orthogonal polynomials 
where the zeros are always real, in the present case the zeros are complex. Szego 0] already 
notes that the zeros of L~' m '(z) are in general non-real except for at most one root. We 
found numerically that all the complex roots live on a smooth curve in the complex plane, we 
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now proceed to obtain the distribution of the roots in the complex plane in the asymptotic 
limit n — > oo. 

Let us define the following Green's function (also known as the Stieltjes transform in the 
mathematics literature): 

An important aspect of our analysis is in the way in which we scale the various arguments 
by n, and normalize by n as above. The remaining variables, such as G, z and y below are 
of 0(1) as n — > oo, whereby various simplifications arise for large n. 
Using the results described in App. (|AD we get: 

„t \ 1 dF(z) 1 d , .. , , 

where, in the present case, we have The usual integral equation 

representation of the associated Laguerre polynomial gives 

1 r exp(f^) 1 fe -n<j>(v,z) 

F(z) = &> — 77 — — dy = & dy, where 

w 2vri 7 (1 - y )-n/d y n+i y 2m J y 

(j) (y, z )=)n(y)-±\n(l-y)-^-. (17) 

The contour can be taken as any closed loop around the origin which does not cross the 
branch line which we take as the real line from y = 1 to oo. In the large n limit the 
integral can be evaluated by a saddle point method, as is standard in statistical physics. 
The saddle point of is determined for each z We can choose the contour to pass through 
the appropriate saddle point. In the usual case, this deformation leads to the asymptotically 
exact result for the "free energy" , here we show that it leads to the exact Green's function. 
The saddle-points are determined through d<f)(y, z)/dy — which gives 



1 — 2d — ndz ± nd J {z — z+) (z — zJ) 

y± = ; — where (18) 

y± 2(1 -d) v ; 

z 1 - 2d i2 [l ~ 
± ird 7i V d 
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Out of the two saddlepoints we choose the one which gives a smaller value for Re[<f)(y, z)] since 
that gives the dominant contribution to the integral. Further it needs to be ensured that it 
is possible to actually draw a contour, enclosing the origin and not crossing the branch-cut, 
such that, along the contour, <p(y, z) takes its minimum value at the saddlepoint. We find 
that either of the branches can be chosen depending on the location of z in the complex 
plane. The condition 

Re[5(j>] = Re[<P(y„, z) - <p(y + , z)] = (19) 

determines the curve in Lz plane on which the roots lie. This condition actually determines 
two curves C\ and C2, where C\ concaves towards the negative Re[z] axis ( like x = —y 2 ) 
and C2 concaves in the other direction. The curves C\ and C2 together form a closed region 
R in Lz plane surrounding the point (1 — 2d)/(-jrd) [see Fig. ([!])]. To the right of the curve 
G\ y the saddlepoint y + dominates. As we cross into the region R, the saddlepoint y_ takes 
over. However when we again cross the curve C2, the saddlepoint y_ continues to determine 
the free energy even though <j)(y + ,z) is smaller. This is because it is no longer possible to 
construct the necessary contour through y + [this is illustrated in Fig. (Q)] and we remain 
stuck to ?/_. Thus it is the curve C\ which determines the locus of roots of the polynomials. 

A brief discussion of the curves is in order. In terms of k — J we replace z by a 
variable ip defined through z = — — k 2 ) + 2/tsinh( , 0)] and hence the branch line of the 
square roots is a curve joining the two points ip = ±^vr/2. The difference in "free energies" 
can be written as 

2 , . ,. 1 + k, 2 , , 1 + Kexp(ib) , 

<ty = z7r + 2V + -cosh(V) h—^g{- T^\}- 20 

k k 1 — Kexp[— ip) 

A particularly nice case is that of d = 1/2 where we have k — 1 and hence another trans- 
formation to b = (b re + ib im ) defined through sinh(6) sinh( , 0) = — 1 gives 

<ty = i7r + 2[&-coth(&)]. (21) 

The roots thus live on the curve b re = Recoth(b re + ibi m ), which simplifies to the curve 
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bim = - arccos [cosh (2 6 re ) — sinh(26 re )/6 re ]. (22) 

Note that z = — 2/(tt sinh(6)), and b re varies between hi = .00079 and hu = 1.19968 the 
values at which the argument of the arccos hits ±1. By choosing the appropriate branch of 
the arccos function we obtain a smooth curve as depicted in Fig. (]T|). 

In Fig. (^) we show curves at different densities, obtained by solving Eq. (|B3) numerically. 
We also plot the roots as obtained from an exact numerical solution of Eq. (JiJ) for n = 54. 



The asymptotic Green's function follows from Eq. (|16l) : 



G{z) = _*Kv>*) 



y± i - y± 



ny± (23) 



dz 

where either of the branches y± is chosen depending on the value of z. A second method of 
deriving the Green's function is described in App. (|B|). 

Density of zeros: The Green's function contains complete information on the roots 
and we now proceed to extract the density of zeros from it. We note that for the case when 
the roots are all located on the real axis, the density of zeros is given by the discontinuity in 
the imaginary part of the Green's function across the branch line formed by the roots. Here 
we develop a generalization of this procedure for the case when the roots are distributed on 
a curve. 

We first assume that the curve can be parametrized by a continuous variable s so that 
points on it are given by z(s) = (x(s),y(s)). The greens function is given by 

n^z - z{si) 
We now define the following limiting function: 

GAs) = lim G(z(s) + rj) = lim - V — ; n = ri x + iri v (24) 

Expanding around si we get: 

lim = lim where fi = — — — 

v^o z(s) - z(si) + Tj v^o z'(si)[(s - Si) +fj] ' z'(si) 

= P- i n 6 -^^ Si gn(f ]y ), (25) 
z'{si) 



where P denotes the principal part having a continuous variation across the curve. Putting 
this into Eq. ( |24| ) and using the definition p(s) = - S(s — s;) we get 

G L {s) = P - m^-sign(fi y ) 
z ys) 

We note that the above result has a simple geometric meaning: the discontinuity is related 
to the variation along the normal to the curve. The discontinuity in the Green's function 
across the curve is thus given by 

Gt(s) - G L (s) = ± 2 J^M = 2n(a + i0), (27) 

z'{s) 

where a and (3 define the real and imaginary parts of the jump. Defining the projected 
densities p x (x) = ^T,S(x - x t ) = p(s)/\x'(s)\ and p y (y) = ~E<% - Vi) = p(s)/\y'(s)\ we 
then get: 

a 2 + (3 2 a 2 + (3 2 

p*( x ) = iffi ; Pviv) = ii ■ ( 28 ) 



In Fig. (^|) we plot the density of zeros, p(s), obtained numerically from Eqns. (|19|) and (p^ 
with the parameter s chosen as the Euclidean length along the curve. 



V. CONCLUSIONS 

In summary we have studied a set of coupled nonlinear algebraic equations which are 
essentially the Bethe ansatz equations for the lowest energy states of the Heisenberg fer- 
romagnetic chain, at small densities. Both a perturbative and an exact solution of the 
equations were obtained. In the former case we find that due to some very special proper- 
ties of the Calogero matrix, we are able to calculate the perturbation series to all orders. 
The exact solution of the equations were obtained following a method due to Stieltjes and 
these are given by the zeros of generalized associated Laguerre polynomials. These zeros are 
distributed on the complex plane and using Green's functions and saddle-point techniques 
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we have obtained the exact asymptotic distribution. It may be noted that there are other 
examples in physics where one needs to study complex zeros of polynomials as for example 
in the case of the Yang-Lee zeros. However this is probably one of the few examples where 
the distribution of zeros in the complex plane has been computed exactly. 



APPENDIX A: 

We give an account of Stieltjes electrostatic interpretation of zeros of the orthogonal 
polynomials [0,3- The case of the Hermite polynomials is widely known in the physics 
literature and has been widely applied. This does not seem to be the case for the other 
orthogonal polynomials. 

We first outline the derivation of the electrostatic interpretation for the Hermite poly- 
nomial. Consider a set of charges in one-dimensions interacting by a logarithmic repulsive 
force and confined within a harmonic potential well. The potential energy of the system is 
given by: 

1 i 

E = - x i ~ 2 E ln ' Xl ~ Xm I ( Al ) 

l=l,n l^m 

The minimum of the energy is given by 

dE/dxi =xi-J2 " = 0- ( A2 ) 

m^l Xl ~ Xm 

Stieltjes result is that the n zeros of the Hermite polynomial H n (x) satisfy this equation. 
To show this, consider the nth-degree polynomial 

g(x) = J] fr- Xl ), (A3) 

l=l,n 

where X\ are the n roots of Eq. (|A2|) . Taking derivatives with respect to x we get 

= e n (* - 9"(x) =e n (?-xm). (A4) 

Hence 



11 



^ = 2£ — =2* (A5) 

since arj satisfy Eq. (|A~2|). Thus we have g"(xi) — 2xig'(xi) = 0. We further note that f(x) = 
g"(x) — 2xg'(x) is a polynomial of degree n with the same roots as g(x) and so they must 
be identical apart from constant multiplicative factor. This constant factor is determined 
by evaluating the coefficient of x n in f(x) which is — 2n. We thus obtain f(x) = —2ng(x) or 

g"{x) - 2xg\x) + 2ng(x) = (A6) 

which is just the equation for the Hermite polynomials. Thus it has been shown that the 
energy function in Eq. ( |Al|) is minimized by a configuration in which the n particles are 
located at the zeros of the H n (x). 

This result can be generalized by replacing 2x\ in Eq. (|A5[) by the expression —q{x{) /p{x{) 
where q(x) and p(x) are arbitrary polynomials of degree 1 and 2 respectively. Taking q(x) = 
q + qix and p(x) = Po + P\X + p 2 x 2 we repeat the previous arguments to obtain 

p{x)g"{x) + q(x)g'(x) + [— n(n — l)p 2 — nqi}g(x) = 0. (A7) 

Denoting the polynomial solutions of this general equation by g n {x) we have the following 
result: The zeros of the polynomial g n (x) correspond to the solution of the set of equations 

q{Xl) ' £ ^— = 0. (A8) 



2p(x l ) ^xi-x 

This corresponds to a minimization of the following n-particle potential energy: 

^E/^ + ^EMsi-^l- (A9) 



2 i J p(xi) 2 ¥m 

We note that these results are valid for arbitrary complex values of the coefficients pi and 
though a physical interpretation cannot be given in all cases. In fact, as we shall now see, 
all the cases which give rise to real confining potentials are precisely the ones corresponding 
to the classical orthogonal polynomials. Physically we expect that for a confining potential 
there will be real solutions to the minimization problem and this is reflected in the fact that 

12 



all orthogonal polynomials of the nth degree have exactly n real zeros (this of-course also 
follows from the Sturm-Liouville theory). We now list the various cases which result in the 
standard orthogonal polynomials. 

(i) V2 7^ 0: In this case, it can be shown that upto linear and scale transformations of 
the variable x, the most general form of the differential equation which leads to a confining 
potential is the following: 

(1 - x 2 )g"(x) + [b - a - (a + b + 2)x\g\x) + n(n + a + b + l)g(x) = 0; a, b > -1 (A10) 

The resulting potential V(x) = — (1+a) In (1 — x) — (l + b) In (1 + x) confines particles within 
the domain — 1 < x < 1. This is the equation for Jacobi polynomials, P^ ,b \x), of which the 
Legendre and Chebyshev polynomials are special cases 0. 

(ii) p2 = 0, pi 7^ 0: In this case the most general form with a bounding potential is 

xg"(x) + (a+ 1 - x)g'(x) +ng(x) = 0; a > -1 (All) 

The resulting potential V(x) — x — (a + 1) lnx confines particles to (0, oo) and corresponds 
to the associated Laguerre polynomials, L®(x) J7J. 

(iii) Finally we have the case pi = P2 = 0, po ^ and this gives rise to Hermite's 
equation which has the harmonic potential V(x) = x 2 /2. 

We note that apart from these special cases the other cases, where Pi, qi are allowed 
to take arbitrary complex values, also may give rise to important physical applications as 
indeed is so for the example of Bethe roots considered in this paper. 

APPENDIX B: 

Here we derive a differential equation satisfied by the Green's function. We also show how 
this also can be applied to the case of the other orthogonal polynomials inorder to extract 
the density of zeros in those cases. We start with Eq. ([UJ) and take a single derivative to 
get: 
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dG(z) 1 .dF.r, 1 d 2 F ,„ . 

"^r^^+Fi? (B1) 

Now F(z) = L~ n l d ~ 1 ( —mrz) and so satisfies the following second order equation (which is 
Eq. (|A11 ) with the substitution y = —mrz and a + 1 = —n/d): 

d? F .dF „ n 
2;——- + (—n/d + imz) — 7m i< =0 

Using this and Eq. ([16]) we can eliminate the derivatives on the right hand side of Eq. ( |B1|) 
to get our equation for the Green's function: 

IdGW _1 £ 

n dz dz z 

In the limit n — > 00, the left hand side vanishes and the roots of the quadratic on the right 
hand give us the two branches of the Green's function. To choose between the branches 
ofcourse requires a examination of the saddle points as described in section (IV). 

This procedure of obtaining the asymptotic Green's function is easily generalizable to 
other polynomial equations and we briefly discuss a few applications to cases where the zeros 
are real. 

(1) Associated Laguerre polynomials with a + 1 = pn > 0: The replacement d = — 1/p 
at once gives us: 



— TXZ — p ± ttJ (z — x+) (z — X-) 
G(z) = — VV A ; 



2z 



with x± = (-p - 2 ± 2^p + 1) /tt. (B3) 

The roots are located in the region x_ < x < x + and the density of zeros, given by the 
discontinuity in the imaginary part of the Green's function is obtained as: 

P(x) = -x){x-xJ) (B4) 

We can also consider the case where the coefficient a does not scale with n but is a finite 
constant. The Green's function and density of zeros are found to be independent of a: 
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2 1 V ' irz s 

p{x) = IJ-L-l; - 4/tt < x < (B5) 
2 y 7r|x| 

(2) Hermite polynomials: In this case we use scaled variables such that the zeros z\ satisfy 
H n (y/nzi) = 0. Proceeding as before we get the following equation for the Green's function: 

ldGjz)_ = _ Q2 + _ 2 

n dz 

The asymptotic Green's function and the density of zeros follow immediately: 



G{z) =z- Vz^2 

p(x) = -V2-x 2 (B6) 

7T 

(3) Jacobi polynomials: If zi satisfy P^ a ' b '(zi) = 0, then the corresponding equation for 
the Green's function is 

1 dG(z) _ 2( \ b — a — (a + b + 2)z . , n + a + b + 1 
n~dz~ = {Z) n(l - z 2 ) W n(l-z 2 ) ' 

We thus get the following asymptotic Green's function and density of states. 

1 



G(z) 



Vz^l 



p{x) = i = ; - K x < 1 (B7) 



7TV 1 — X 2 



The results in Eq. (P5|) , (p^ ) and ( |B7| ) are identical to those given in Refn. |J where they 



have been obtained by different methods. 
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FIGURES 

FIG. 1. The curves C\ and C2, in the complex z plane, for density d = 1/2. At this density 
the curves are explictly given by the formula in Eq. (22). The curve C\ corresponds to the true 
line of singularities on which the roots lie. There is also a change of branch from y + to y~ as we 
cross the vertical lines C% from right to left. But there is no discontinuity of the Green's function 
across these lines and hence there are no zeros on these lines. 
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FIG. 2. The figure shows, schematically, contours of constant </>+ and in the complex y 
plane (with d = 0.6 and z = —0.4 + 0.3i). For this case 0+ < 0-. It is clear that one cannot draw 
the appropriate contour through the point y + and so we are forced to choose the point y_ even 
though it has a higher value of 4>. This illustrates why there is no change in branch when we cross 
the curve C2. 
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FIG. 3. The solid curves give the asymptotic distribution of zeros in the complex z-plane 
obtained from a numerical solution of Eq. ( |i~9|) for d = 0.2,0.4,0.6 and 0.8 ( extreme left ). Also 
shown are numerically obtained roots (points) for n = 54. 
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FIG. 4. The asymptotic density of zeros p(s) as a function of the Euclidean length along the 
curve for d = 0.2, 0.4, 0.6 and 0.8 . The densities are evaluated numerically using the formula in 
Eq.(||). 
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